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Abstract 

Background: Fasting induces specific molecular and metabolic adaptions in most organisms. In biomedical 
research fasting is used in metabolic studies to synchronize nutritional states of study subjects. Because there is a 
lack of standardization for this procedure, we need a deeper understanding of the dynamics and the molecular 
mechanisms in fasting. 

Results: We investigated the dynamic changes of liver gene expression and serum parameters of mice at several 
time points during a 48 hour fasting experiment and then focused on the global gene expression changes in 
epididymal white adipose tissue (WAT) as well as on pathways common to WAT, liver, and skeletal muscle. This 
approach produced several intriguing insights: (i) rather than a sequential activation of biochemical pathways in 
fasted liver, as current knowledge dictates, our data indicates a concerted parallel response; (ii) this first 
characterization of the transcriptome signature of WAT of fasted mice reveals a remarkable activation of 
components of the transcription apparatus; (iii) most importantly, our bioinformatic analyses indicate p53 as central 
node in the regulation of fasting in major metabolic tissues; and (iv) forced expression of Ddit4, a fasting-regulated 
p53 target gene, is sufficient to augment lipolysis in cultured adipocytes. 

Conclusions: In summary, this combination of focused and global profiling approaches provides a comprehensive 
molecular characterization of the processes operating during fasting in mice and suggests a role for p53, and its 
downstream target Ddit4, as novel components in the transcriptional response to food deprivation. 
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Background organs switch from glucose to fatty acids as the primary 

Fasting processes are evolutionarily highly conserved energy source. These fatty acids are mainly supplied from 

adaptive responses to food deprivation in all organisms. adipose tissue stores where they are released from trigly- 

Mammals follow a day/ night pattern that includes a ceride droplets by the process of lipolysis [2]. However, 

daily resting period during which the body is in a fasted the brain is not equipped to derive energy from fatty acids, 

state and must adapt to the lack of nutrient intake by but uninterrupted maintenance of its function is essential 

changing the metabolic state in several organs and at the to organismal survival. Hence, the perhaps most astound- 

systemic level. When deprived of dietary nutrients, the ing metabolic change in fasting is the massive increase in 

body initially derives glucose from glycogen stores, which the hepatic production of ketone bodies, which the brain 

are quickly depleted [1]. If fasting continues, peripheral can use as a secondary energy source [3]. The emergence 

of biochemical pathways that produce and then utilize ke- 
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While this is crucial for survival of millions in poor, 
under-developed countries, Western societies are facing a 
pandemic of obesity with more than one third of the 
population being obese in countries like the United States 
of America [4]. In this context, fasting is still one of the 
most efficient weight loss measures and therefore a 
valuable therapeutic tool for the management of obesity 
[5]. Indeed, caloric restriction and alternate-day fasting 
have both been shown to increase longevity [6] and reduce 
metabolic risk factors [7]. Further, a ketogenic state as in 
fasting can also be elicited by low-carbohydrate diets 
which are among the most efficient weight loss programs 
[8] and are also used in clinical settings to treat forms of 
epilepsy [9]. For these reasons it is imperative to add to 
the current knowledge about the molecular underpinnings 
and the systemic consequences of fasting. 

In experimental animals, fasting is frequently used in 
studies where nutritional states need to be synchronized 
in order to keep biological variation (introduced by vary- 
ing patterns of food intake) minimal and, thereby, to pro- 
vide a steady baseline for all measured parameters. The 
major problem with this approach is that there is no 
standard protocol that the research community adheres to 
and studies report fasting regimes in a wide range (from a 
few hours to days) and with different day-time starting 
points, often disregarding the intricate circadian regulation 
on the organ and the systemic level [10]. Consequently, a 
systematic meta-study identified fasting as one of the main 
sources of variation between different transcriptomic 
studies [11]. Hence, to estimate the impact on measured 
experimental parameters when animals are fasted, a 
deeper understanding of this process is needed. For in- 
stance, whereas the fasting response of several mouse 
tissues (liver, brain, gut, muscle, kidney) has been inves- 
tigated at the transcriptome [12,13] and proteome [14] 
level, a comprehensive view on gene regulation in white 
adipose tissue of fasted mice is still missing. 

In this work we investigated the dynamic and circadian 
responses to a fasting stimulus by measuring serum 
parameters and liver gene expression in fasted mice at 



several time points. Moreover, we measured the global 
transcriptome response to fasting in white adipose tissue, 
liver, and skeletal muscle utilizing microarray technology 
followed by a host of bioinformatic analyses. Interesting 
outcomes from the fasted adipose tissue data include a 
strong enrichment of upregulated genes coding for pro- 
teins that regulate transcription as well as apoptosis- 
related genes. A major finding is the identification of the 
p53 signaling pathway as a common mediator of the fast- 
ing response in all three investigated tissues. Driven by 
these results we focused on the p53-target gene Ddit4, 
which was upregulated by fasting in all three tissues. Our 
experiments showed that, in cultured adipocytes, Ddit4 is 
inducible by p53 activation and its ectopic expression is 
sufficient to augment lipolysis. Therefore we describe a 
new molecular component in the fasting-response down- 
stream of p53. 

Results and discussion 

Kinetics of serum parameters and liver gene expression in 
mice during a 48 hour fasting period 

Despite a considerable amount of knowledge about the 
mechanisms of fasting response in different organisms 
(e.g. [1,12]), its exact timing, especially at time points 
earlier than 12 hours, is still unclear. To assess the timely 
fasting response of serum metabolites we measured blood 
glucose, non-esterified fatty acids (NEFA), glycerol, and 
(3-hydroxybutyrate, as well as the hormones insulin and 
corticosterone, at the start of the experiment and at five 
time points during a 48 hour period (0, 3, 6, 12, 24 and 
48 hours) both in fasted and in control (i.e. ad-libitum 
fed) male C57B1/6 J mice. The experiment was started 
at the beginning of the light cycle (9 a.m.) when mice 
are in their inactive phase. For these time points we 
determined body weight and weight changes within 
each group (Table 1) and measured the expression of 
genes in pathways central to the fasting response in liver, 
which is considered the main organ for maintaining 
systemic energy homeostasis. 



Table 1 Weight changes of study animals 




Absolute weight [g] 


Weight change to initial weight within groups [g] 


ZT 


Fed (n = 5) 


Fasted (n = 5) 


Fed (n = 5) 


Fasted (n = 5) 


2 Oh 


26.7 ± 0.6 




-0.21 ±0.62 




5 3 h 


28 ± 0.74 


26.2 ± 0.65 


-0.1 3 ±0.1 4 


-2.1 3 ±0.2 


8 6 h 


25.5 ± 1.73 


28 ± 0.68 


-1.14±0.11 


-2.26 ±0.1 9 


14 12 h 


26.8 ± 0.93 


24.8 ±0.57 


-1.09 ±0.15 


-2.96 ±0.1 4** 


26 24 h 


25.1 ±0.69 


22.9 ± 0.33 


-0.48 ±0.1 5 


-3.5 ±0.14** 


50 48 h 


26.9 ± 0.44 


22 ± 0.55*** 


-0.66 ± 0.29 


-5.9 ± 0.42*** 



Significant differences were determined with a 2-way ANOVA followed by a Bonferroni posttest to determine significance for each time point (in bold: **p < 0.01; 
***p < 0.001). ZT: zeitgeber time (hours). 
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Immediate upregulation of hepatic gluconeogenesis by 
fasting 

A well-known and essential transcriptional response to 
fasting is the upregulation of Pckl and G6pc mRNAs 
[15-17]. Pckl (all gene names referring to mRNA tran- 
scripts in this article are in italic and given according to 
the unique NCBI gene symbol nomenclature; full names 
can be found in the abbreviations list) participates in 
controlling the flux into the gluconeogenesis (GNG) 
pathway, while G6pc catalyzes the last step in this path- 
way converting glucose-6-phosphate to glucose. Activa- 
tion of liver GNG is essential in times of undernutrition 



to provide the brain with glucose, its primary energy 
source. We found that this upregulation occurs already 
after 3 hours fasting and continues, at least in case of 
Pckl mRNA, throughout 48 hours of fasting (Figure 1A 
top row). Further, we measured continuous upregulation 
of Pcx as well as Gyk (Figure 1A top row). Pcx converts 
pyruvate to oxaloacetate, which in turn serves as substrate 
for Pckl. Gyk, which has already been shown to be upreg- 
ulated by fasting in a Ppara-dependent manner [18], cata- 
lyzes the first step in the glycerol phosphate shuttle which 
utilizes glycerol released from triglyceride stores to be 
converted into lipids or shunted into the GNG pathway. 
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Figure 1 Kinetics of selected metabolites/hormones and expression of liver genes during a 48 hour fasting period. (A) and (B) Serum 
parameters and liver mRNA levels were measured in 10-12 weeks old fasted and ad-libitum fed male C57BI/6 J mice at several time points (each 
n = 5). Shading indicates light or dark cycle. Relative expression values from qPCR measurements were normalized to 36b4 expression and the 
values for time point 0 were set to 1. Significance was determined with a 2-way ANOVA followed by a Bonferroni posttest to determine 
significance for the single time points (*p < 0.05; **p < 0.01; ***p < 0.001). 
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In accordance with the early upregulation of gluconeo- 
genic genes, serum levels of corticosterone (glucocorti- 
coids are known as positive regulators of hepatic GNG by 
regulating Pckl and G6pc levels via the glucocorticoid 
receptor [19,20]) rose immediately after onset of fasting, 
showing >4-fold higher levels later in the fasting period 
(Figure 1A bottom row). However, during the first 6 hours 
of fasting, blood glucose levels are similar in fasted and 
control-fed mice despite an immediate drop in serum in- 
sulin in fasted mice (Figure 1A bottom row). This might 
be explained by an early activation of glycogenolysis in the 
major glucose storing tissues (i.e. muscle and liver [21]). 
After 12 hours, glycogen stores seem to be depleted and 
serum glucose levels are significantly lower in fasted mice, 
suggesting that liver GNG is compensating only partly for 
the lack of dietary glucose by keeping blood glucose levels 
around 100 mg/dl (Figure 1A bottom row). Hence, our 
data on the expression profiles of gluconeogenic genes 
show an immediate upregulation of hepatic GNG 3 to 
6 hours after food withdrawal, despite a delayed reduction 
of serum glucose levels. 

Increase of liver ketogenesis by three hours after onset of 
fasting 

One substrate for hepatic GNG is glycerol, which is 
mainly derived from adipose tissue lipolysis where triglyc- 
erides are hydrolyzed to glycerol and NEFA. As indicated 
in Figure IB, this process results in an increase in serum 
glycerol and NEFA after 6 hours of fasting. Serum NEFA 
levels normalize to the control-fed levels at the 24 hour 
time point, presumably due to increased uptake and 
utilization of free fatty acids in several tissues (mainly liver 
and skeletal muscle) while glycogen stores are further de- 
pleted. Apart from being a major energy source in times 
of nutrient deprivation, adipose tissue-derived fatty acids 
serve as substrate for ketogenesis, the synthesis of ketone 
bodies (mainly p-hydroxybutyrate and acetoacetate) in the 
liver. These ketone bodies can be used as secondary en- 
ergy source by the brain, which cannot utilize fatty acids 
directly [3]. In our data set, serum (3-hydroxybutyrate 
levels are increasing steadily over 48 hours of fasting com- 
pared to ad libitum-fed control mice with the first signifi- 
cant increase at 3 hours (Figure IB bottom row). This 
increase is concordant with the upregulation of liver 
Hmgcs2 mRNA which is a key enzyme in liver ketogenesis 
[22] (Figure IB top row). Interestingly, liver ketogenesis 
and serum p-hydroxybutyrate are increased hours before 
blood glucose levels begin to drop (Figure 1A bottom 
row). This dynamic argues for the involvement of other 
sensors (for instance gut hormones or insulin and gluca- 
gon signaling) detecting the absence of nutritional carbo- 
hydrates and signaling to liver to upregulate ketogenesis 
early in fasting. A recent report introduced an intriguing, 
novel functional role for ketone bodies. By inhibiting 



histone deacetylases, p-hydroxybutyrate was shown to in- 
duce the expression of genes that protect against oxidative 
damage in a variety of tissues [23] . In that study, a change 
in histone deacetylation in kidney was shown at serum 
p-hydroxybutyrate concentrations greater than 0.5 mM, 
a level we measured already after 3 hours of food with- 
drawal. Thus it is conceivable that this p-hydroxybutyrate- 
mediated mechanism could be activated early on, to set 
the stage for transcriptional regulation in response to 
fasting. 

Early fasting-onset upregulation of Ppara target genes 

Fgf21 is a recently discovered hormone shown to be a 
major hub in the hepatic response to fasting by regulat- 
ing fatty acid and glucose metabolism [24-26]. In our 
data Fgf21 mRNA was elevated in liver of fasted animals 
throughout the 48 hours, peaking at 24 hours (Figure IB 
top row). Furthermore, Fgf21 was shown to mediate its 
effect partly via upregulation of Ppargcla [27], a tran- 
scriptional coactivator we found to be highly increased 
by fasting in accordance with Fgf21 levels (Figure IB 
top row). Ppargcla in turn increases expression of many 
fasting response genes by binding and coactivating 
transcription factors such as Ppara and glucocorticoid 
receptor [15]. Along these lines, most genes shown in 
Figure 1 are Ppara target genes {Pckl, G6pc, Pcx, Gyk, 
Hmgcs2, Fgf21; see [18,25,28-30] and references therein) 
arguing for the central role of this transcription factor 
during fasting, evident from the phenotype of fasted 
Ppara knock-out mice [31]. However, the modest changes 
of liver Ppara mRNA levels (Figure IB top row) are un- 
likely to cause the strong alterations in Ppara targets 
(compare for instance to Fgf2l, a functional Ppara target 
[25]). Rather the transactivation of Ppara by endogenous 
ligands (fatty acids and their derivatives [32]), coactivation 
by Ppargcla, and synergistic regulation by other fasting- 
regulated transcription factors (e.g. Foxa2 [33] and Creb 
[34]) could lead to the magnitude of increase of its target 
genes. 

In summary, comparing expression of key liver fasting 
genes to serum parameters shows a coherent picture 
suggesting, in accordance with other recent studies 
[12,13], a parallel activation of fasting-induced pathways 
rather than a sequential response as historically believed 
[1]. This response is activated as early as 3 to 6 hours 
after food withdrawal and reaches a steady state be- 
tween 12 and 24 hours. Our data further underlines that 
Ppara acts as one major fasting hub, by coordinating 
expression of its target genes. Hence, we provide a 
detailed view of molecular response kinetics during a 
48 hour fasting period in mice allowing one to extrapo- 
late on the timely regulation of the fasting response in 
liver and in the whole organism. 
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Global changes in transcriptome signatures of white 
adipose tissue, liver, and skeletal muscle in fasted mice 

Next, we aimed to elucidate RNA abundance responses 
to fasting in a systematic and genome-wide manner. 
Most of the parameters determined in Figure 1 show the 
highest difference between fasted and fed states at 
24 hours initiating the experiment, suggesting that the 
metabolic adaption to fasting has reached a first steady 
state. For this reason we chose this time point for tran- 
scriptome analyses of epididymal white adipose tissue 



(WAT), liver (LIV), and skeletal muscle (SM) of fasted 
and control mice (Figure 2A). Tissue-derived mRNA 
from five fasted and five control mice were hybridized to 
Affymetrix GeneChip arrays (Mouse Gene 1.1 ST). Two 
hybridizations (fed liver #4 and fasted skeletal muscle 
#1) were outliers as determined by principal component 
analysis and therefore excluded from further analysis. 
Hierarchical clustering of the remaining microarray data 
sets showed that experiments strongly cluster by tissue 
type (Figure 2B). The heatmap contains about 7000 
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Figure 2 Overview of microarray experiments in WAT, LIV, and SM 24 hours after onset of fasting. (A) Experimental design of the 
transcriptome study (WAT = white adipose tissue, LIV = liver, SM = skeletal muscle, PCA = principal component analysis). (B) The heatmap 
contains about 7000 probe sets differentially expressed > 1.3-fold (FDR5) between fasted and fed in at least one condition. Hierarchical clustering 
of experiments clusters biological replicates together. Two outlier replicates were identified by principal component analysis and removed from 
the data set. (C) Numbers of up- and downregulated (1.3x, FDR5) microarray probe sets and RefSeq -annotated genes (in parentheses) in WAT, LIV, 
and SM of mice fasted for 24 hours. Additional file 1 provides detailed expression values for these genes. 
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probe sets differentially expressed by > 1.3-fold at a false 
discovery rate < 5% after multiple testing adjustment 
(FDR5) in at least one condition. WAT showed the highest 
number of differentially expressed probe sets between fast- 
ing and fed states (1850 upregulated and 1712 downregu- 
lated upon fasting, collapsing into 1491 upregulated and 
1469 downregulated RefSeq-annotated genes), more than 
double the number of genes regulated in liver (Figure 2C). 
Additional file 1 provides the lists of probes expressed 
more than 1.3-fold with an FDR5 for each of the three 
tissues. 

Functional annotation clustering of genes regulated in 
fasted white adipose tissue 

Although the transcriptome signatures of other mouse 
organs, including liver and muscle, during fasting were 
extensively investigated [12,13,34] and analyses of adipose 
tissue transcriptomes during fasting are available (rats 
[35,36], chicken [37], pigs [38], lactating goats [39]) the 
current study is, to our knowledge, the first to focus on 
characterization of the global gene expression response of 
adipose tissue in fasted mice. Hence, we focused our initial 
analysis on the WAT dataset and performed functional 
annotation of the genes up - or downregulated in fasted 
WAT (1.3x, FDR5) by mapping these lists onto gene 
ontology (GO; [40]) and KEGG pathways [41] utilizing the 
DAVID tool [42,43] . To visualize the results we employed 
the functional clustering option, which combines redun- 
dant entities, based on similarity of gene lists, into clusters. 
Figure 3A and B respectively show the resulting clusters 
from the lists of "WAT downregulated" or "WAT upregu- 
lated" genes together with contributing GO domains in 
parentheses. The x-axis shows the cluster group en- 
richment scores representing a geometric mean (-log) 
of p-values of entities in each cluster. We only show 
clusters with an enrichment score larger than 3 be- 
cause they contain at least one entity with a significant 
p-value after multiple testing correction (<0.05). The 
detailed results of the functional clustering including 
terms/pathways contributing to each cluster are given 
in Additional file 2. 

Functional clustering of genes downregulated by 
fasting in WAT (Figure 3A) yields terms related to 
"mitochondria" in the cellular component (CC) GO cat- 
egory, indicating suppression of mitochondrial processes 
in times of energy scarcity. In the GO domain molecular 
function (MF) we find "G TP-associated" terms and GO 
terms related to "nucleotide binding". GO biological pro- 
cesses (BP) and KEGG pathways are found in clusters re- 
lated to "steroid and sterol metabolism", "carbohydrate 
metabolism" and "biosynthesis of fatty acids and amino 
acids". Downregulation of steroid and sterol metabolism is 
analyzed and discussed in more detail later. The cluster of 
carbohydrate metabolism is mainly comprised of GO 



terms referring to catabolism of monosaccharides, like 
hexose or glucose, and might reflect the shift from glucose 
usage to fatty acid (3-oxidation in prolonged fasting. Also 
the downregulation of fatty acid and amino acid biosyn- 
thesis is a plausible reaction of a fasting cell or tissue that 
needs to suppress energy consuming processes which 
increase anaplerotic pressure. 

Genes upregulated in fasted white adipose tissue are 
strongly connected to transcriptional regulation 

Genes that are upregulated in fasted WAT (Figure 3B) 
show the highest enrichment score for the GO MF cat- 
egory related to "ion binding", which is a rather general 
and unspecific term. However, we find a cluster con- 
taining 127 genes that relate to "apoptosis" and "cell 
death" indicating that apoptotic pathways are activated 
during fasting in WAT. Finally, we find high enrichment 
in a cluster containing "nucleus"-related GO CC terms 
as well as a cluster with terms referring to "regulation of 
transcription" (BP and MF), both of which indicate a 
strong positive influence of fasting on the transcrip- 
tional program of adipose tissue. This is consistent with 
WAT being the tissue with the highest number of differ- 
entially regulated genes when compared to LIV and SM 
(Figure 2C). To substantiate this finding we tested 
whether the genes in the upregulated list are enriched 
for transcriptional regulator molecules such as tran- 
scription factors, cofactors and chromatin remodelers 
as defined by the manually curated AnimalTFDB [44]. 
As a control we also tested for enrichment in the list of 
genes downregulated by fasting in WAT which did not 
map to "regulation of transcription"-associated GO terms 
in our DAVID analysis. As shown in Additional file 3 the 
WAT upregulated list contains at least double the number 
of transcriptional regulators, compared to the WAT 
downregulated list (e.g. 122 vs. 50 transcription factors), 
and comprises 9%, 12%, and 16% of annotated mouse 
transcription factors, cofactors, and chromatin remode- 
lers, respectively. Additional file 3 also lists the expression 
values and the annotations of regulated genes in WAT 
according to the transcription factor subclasses defined in 
AnimalTFDB. Based on ranking of the entire data set, the 
gene-set enrichment analysis tool [45] performs a similar 
analysis and yields a significant enrichment for the sub- 
categories (all but HOX factors) of transcription factors 
from AnimalTFDB for upregulated genes (see Figure 3C 
and Additional file 3). These transcription factor families 
contained 43 genes encoding zinc finger transcription 
factors (including kruppel-like factors Klf4, Klf9, and 
Klfl5), 9 genes encoding basic leucin zippers (bZIP), 
and 6 forkhead box genes that are upregulated in fasted 
WAT (see Additional file 3). Among the transcription 
factor-encoding genes with the highest upregulation by 
fasting (14-fold) we find Irf4, which has been shown to 
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Figure 3 Functional annotation of genes regulated in WAT 24 hours after onset of fasting reveals activation of transcription regulatory 
components. (A) and (B) DAVID functional clustering was performed with genes upregulated (A) or downregulated (B) by fasting (1.3x, FDR5). 
Shown are the numbers of unique genes (inside of bars) and the enriched clusters (left of bars with contributing DAVID domains in parentheses). 
The x-axis shows the cluster enrichment scores representing a geometric mean (-log) of p-values of entities in each cluster. Clusters with an 
enrichment score larger than three are shown. CC = cellular components, MF = molecular functions, BP = biological processes. (C) GSEA enrichment 
plots are generated by ranking all genes in the dataset by expression values, mapping transcription factor (TF) sub-categories (bZIP = basic 
leucine zippers; Fox = forkhead box), cofactors, and chromatin remodelers, and calculating a cumulative enrichment score (green line). 
NES = normalized enrichment score, FDR = false discovery rate. 



tip the scales between lipogenesis and lipolysis in the 
latter direction in a fasting-induced manner [46]. Con- 
sequently, fat-specific Irf4 knock-out mice are deficient 
in lipolysis and show increased adiposity [46]. Further- 
more, genes encoding the transcription factors Ziml 
and Peg3 are massively upregulated by fasting (21 -fold 
and 6-fold, respectively). Interestingly, these two genes 
are adjacently located at an imprinted region on mouse 
chromosome 7. The Peg3 knock-out mouse model de- 
velops increased adiposity despite lower food intake. 
This was attributed to developmental deficiencies that 
lead to aberrant leptin signaling in the hypothalamus 
[47]. However, our data suggest a direct involvement of 
the Peg3/Ziml locus in adipose tissue biology. 



The present study is, to our knowledge, the first to 
focus on the characterization of the transcriptome re- 
sponse to fasting of WAT in mice. Surprisingly, we do 
not find GO mapping to the biological process "lipid 
catabolism" which would be expected to be prominent 
during fasting in WAT. The absence thereof could be 
explained by abundant posttranscriptional regulation 
of lipolysis by kinases such as protein kinase A [48] 
and AMP -activated protein kinase [6] which are not 
reflected at the transcript level. Instead our analyses 
specifically reveal an unexpected upregulation of cell 
death pathways as well as a strong enrichment of tran- 
scriptional regulators among genes activated by fasting 
in WAT. 
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p53 signaling as top ranking pathway in the fasting 
response of major metabolic tissues 

By producing and exporting glucose, fatty acids, gly- 
cerol, and ketone bodies WAT, LIV, and SM represent 



the organs mainly responsible for energy homeostasis 
during a fasting period in mammals [1,21]. To reveal 
common pathways we focused on 200 genes that were 
regulated by fasting in all three tissues (Figure 4A, 
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Figure 4 Analyses of differentially expressed genes regulated by fasting in all three tissues revealing an involvement of p53 and Srebp 
signaling pathway. (A) Venn diagram of genes differentially regulated (1.3x, FDR5) in WAT, LIV, and SM and their intersections. (B) Two hundred 
genes differentially regulated in all three tissues were submitted to DAVID functional annotation using GO terms and KEGG pathways. Count 
shows the number of mapped genes and Benjamini-Hochberg's correction was used for adjusting p-value for multiple testing. (C) Metacore 
mapping as obtained from analysis with the commonly regulated gene list. The thermometers display expression levels of mapped genes 
(1: WAT, 2: SM, 3: LIV). (D) Heatmap of p53 target genes discussed in the text. Values are log2 transformed. White dots mark highest expression 
between the 3 tissues. (E) and (F) Manually compiled heatmaps indicating expression levels of genes comprising the FA (E) and cholesterol 
biosynthesis (F) pathways. Values are log2 transformed. White dots mark highest expression between the 3 tissues. 



Schupp et al. BMC Genomics 2013, 14:758 
http://www.biomedcentral.com/1471 -21 64/1 4/758 



Page 9 of 20 



termed "common list" in the following). Table 2 lists 
these genes ranked by decreasing average expression 
level Mapping the common list to GO biological pro- 
cesses and KEGG using DAVID yielded one KEGG 
pathway as significantly enriched after multiple testing 
correction: the p53 signaling pathway (Figure 4B). Add- 
itionally, an independent analysis with Metacore (GeneGo, 
Inc.) focusing on networks overrepresented in the com- 
mon list reveals the p53 node as the second highest scor- 
ing network hub (see Additional file 4 for genes that 
constitute this network; the highest scoring hub was the 
ubiquitous transcription factor Spl). The p53 transcription 
factor is known to regulate a number of tumor suppressor 
pathways including cellular senescence, apoptosis, and 
DNA repair in response to various stressors [49,50]. About 
50% of human cancers carry mutations in the p53 gene 
[51]. However, evidence is now accumulating that p53 
plays a prominent role in (lipid) metabolism as well 
[52,53]. There is no direct evidence of a physiological 
importance of p53 signaling in fasting, but some studies 
report an activation of p53, and an induction of its target 
genes, upon glucose deprivation in cultured cells [54,55]. 
In our dataset we observe a host of p53 target genes ap- 
parently regulated by fasting, and a heatmap of the genes 
described below is shown in Figure 4D. For instance, the 
canonical p53 target Cdknla [49], a cell cycle regulator 
better known as p21, is one of the top ranked genes in 
Table 2 and is the gene with the highest fasting- induced 
upregulation in liver (>10x). Another example is lipin-1 
encoded by the gene Lpinl, which has recently been 
described as a novel p53 target [55]. In our data Lpinl is 
upregulated in WAT, LIV, and SM by fasting. Depending 
on its subcellular location lipin-1 can act to enhance fatty 
acid oxidation by interacting with Ppara and Ppargcla (in 
the nucleus) [56] or, as a phosphatidate phosphatase in the 
cytosol, to perform a key step in triglyceride biosynthesis 
[57]. Being upregulated by fasting suggests activity of its 
nuclear form to drive oxidation of fatty acids in lipid 
storing tissues. Other p53 targets are the Sestrins, Sesnl 
(upregulated by fasting in WAT and SM) and Sesn2 
(upregulated in all 3 tissues), which have recently been 
shown to be induced by fasting in liver and to protect it 
from oxidative damage in a fasting/refeeding scenario 
[58]. Interestingly, Sen et al. report that p53 physically 
interacts with Ppargcla [59], which we show to be 
strongly upregulated in liver (Figure IB) by qPCR as 
well as in our microarray study in WAT and LIV (it is 
downregulated, however, in SM). In their study Sen 
et al. show that Ppargcla can bind and thereby direct 
p53 to the promoters/enhancers of pro-arrest as well as 
metabolic target genes. Thus, this interaction of p53 
and its coactivator Ppargcla could lead to a tissue- 
specific coordination of p53 to target genes relevant for 
a proper response. 



Fasting-mediated downregulation of Srebp pathways is 
common to major metabolic tissues 

It was reported that p53 is capable of suppressing the 
promoter of the gene encoding the transcription factor 
sterol-regulatory element binding protein 1 (Srebpl 
[60]). The Srebp family consists of three members: 
Srebpla and Srebplc, both transcribed from the same 
gene Srebfl, with Srebp- lc being the predominant 
isoform to regulate lipogenesis in metabolic tissues, 
such as WAT and LIV; and Srebp2, transcribed from 
the Srebf2 locus and responsible for regulation of sterol 
metabolism [61-63]. Consistent with these reports, a 
Metacore analysis on the common list delivers "Regulation 
of lipid metabolism" as top ranking pathways with the 
Srebp family of transcription factors in its center 
(Additional file 5). As mapped in Figure 4C, we find 
Srebfl and Srebf2 downregulated in all three tissues 
(LIV Srebp2 is not mapped because it shows a significant 
1.26x downregulation, not meeting our 1.3x cut-off). Con- 
cordantly, established Srebpl downstream genes, coding 
for enzymes which encompass the fatty acid biosynthesis 
pathway like Acss2, Acaca, Fasn, Scdl, and Scd2 [64], are 
robustly down-regulated by fasting in our data (Figure 4E, 
Table 2). Further, DAVID analyses for the set of commonly 
regulated genes (Figure 4B) as well as for regulated liver 
genes (Additional file 6) yields GO biological processes 
that refer to steroid metabolism and, more specifically, to 
cholesterol biosynthesis. As a pivotal regulator of choles- 
terol homeostasis in cells, Srepb2 mediates its effects by 
control of de novo synthesis and/or by regulation of chol- 
esterol import [65]. To regulate de novo synthesis, Srebp2 
(and to some extent also Srebpl) binds to promoters of 
most of the enzymes in the cholesterol biosynthesis 
pathway [63]. We find that transcripts for most of these 
enzymes are significantly downregulated or show at 
least a trend to downregulation (Figure 4F). The intern- 
alization of cholesterol is regulated by Srebp-mediated 
transcription of the LDL-receptor Ldlr (down in all 
three tissues). Figure 4F shows a heatmap of these 
Srebp2 target genes [63] and impressively demonstrates 
reduction of transcripts for nearly all components of the 
cholesterol biosynthesis and uptake pathway in the in- 
vestigated tissues. Given the extensive and intricate net- 
works that, depending on sterol or nutrient availability, 
post-translationally regulate processing and activation 
of membrane-bound SREBPs [66], our finding that fast- 
ing broadly regulates Srebp-dependent pathways already 
at the transcriptional level is rather surprising. Espe- 
cially in liver, this downregulation of cholesterol biosyn- 
thesis in combination with the upregulation of Hgmcs2 
(Figure 1A and Figure 4F), which condenses acetoacetyl 
CoA and acetyl CoA to (3-hydroxy (3-methylglutaryl CoA, 
hints to a diversion of substrate from sterol synthesis to 
ketone body synthesis. 
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Table 2 Fold difference in expression (fasted/fed) of transcripts regulated in WAT, LIV, and SM by fasting 
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Table 2 Fold difference in expression (fasted/fed) of transcripts regulated in WAT, LIV, and SM by fasting (Continued) 
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Table 2 Fold difference in expression (fasted/fed) of transcripts regulated in WAT, LIV, and SM by fasting (Continued) 
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Genes with >1.3x change (FDR5) in all investigated tissues. 



Hence, based on our analyses on genes regulated in 
WAT, LIV, and SM, we hypothesize a shared mechanism 
that responds to the fasting stimulus in all three tissues: 
The p53 signaling pathway is activated by fasting in WAT, 
LIV, and SM (perhaps via AMPK as described in [67]). 
Coactivators like Ppargcla direct p53 to promoters/ 
enhancers of genes targeted for transcription activation 
or repression when nutrients are lacking. Upregulation 
of p53 targets such as Lpinl might contribute to the 
shift of fasted tissues to fatty acid oxidation to provide 
energy substrates. The observation that p53 knock-out 
mice are incapable of inducing liver fatty acid oxidation 
upon fasting [54] underlines this hypothesis. In addition, a 
p53-mediated downregulation of Srebpl is followed by a 
repression of fatty acid biosynthesis (Figure 4C and E). 
However, we note that other (parallel) pathways that are 
regulated at the post-translational level, and therefore not 
reflected in RNA levels, could be responsible for some of 
the observed effects, such as downregulation of Srebp 
transcripts. One conceivable example is the activation of 
AMPK by fasting which, via subsequent deactivation of 
mechanistic target of rapamycin complex 1 (mTORCl), 
could be responsible for the decrease in Srebp mRNA 
[62]. However, based on our analyses we propose a novel 
and potentially crucial role for p53 in fasting, which even- 
tually could manifest in profound transcriptional changes 
in several metabolic pathways. Although functional proof 
of this mechanism is necessary, several reports support 
our model [52,54,55,58-60]. 

Validation of expression of top-ranked genes commonly 
regulated by fasting in WAT, LIV, and SM 

To validate the microarray data by means of qPCR we se- 
lected the top 3 genes from Table 2 (highest average ex- 
pression across WAT, LIV, and SM: Ddit4, Cdknla, and 
Cebpd), none of which have previously been functionally 
described in the context of fasting in mice, as well as Perl 
and Fasn, known responders to food deprivation and 



thereby positive controls (Figure 5). Perl is an essential 
regulator in the core clock machinery of circadian rhythm 
[10,68] and, in LIV and AT at the beginning of the light 
phase, it was previously shown that its expression is higher 
in fasting compared to control fed mice [69,70]. This dif- 
ference is due to a circadian phase shift that takes place 
during an extended fasting period [69]. As mentioned 
above Fasn has long been known to be a downstream tar- 
get gene of Srebpl and is downregulated due to the drop 
in Srebpl levels during fasting [71,72]. The basic-leucin 
zipper transcription factor Cebpd has been described in a 
number of cellular contexts, such as osteogenesis and adi- 
pogenesis (reviewed in [73]). Cebpd expression is known 
to respond to glucocorticoids and to increased cAMP 
levels, both of which could explain its upregulation upon 
fasting. Cdknla (p21) as a major p53 target gene, is mainly 
described as a cell cycle and apoptosis regulator that in- 
hibits cyclin-dependent kinases [74] and has no known 
role in fasting. Finally, Ddit4, a gene initially reported to 
be readily induced by dexamethasone [75] as well as upon 
certain cellular stresses (e.g. hypoxia [76] and exercise in 
muscle [77]), shows the highest extent of upregulation in 
WAT and SM in the common list (Table 2). Interestingly, 
it has also been described as a p53 target gene [78], which 
led us to further investigate it (see below). Hence, we con- 
firmed fasting-mediated regulation of all genes selected 
for qPCR validation in all three tissues (Figure 5A) and 
show a strong correlation with the microarray measure- 
ments for all these genes (Figure 5B, r 2 = 0.91). This intro- 
duces three intriguing and novel players in the response to 
fasting. 

Ddit4 is fasting-induced in WAT, LIV, and SM and is 
inducible by p53 activation in cultured adipocytes 

To investigate p53 signaling as a common fasting regula- 
tor in WAT, LIV, and SM, we focused on DNA damage- 
induced transcript 4 (Ddit4, aliases: Reddl, Dig2, Rtp801), 
the top ranking gene in Table 2 which has been described 
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Figure 5 Validation of selected candidates by qPCR. (A) Five 
genes differentially regulated in WAT, LIV, and SM were selected for 
qPCR validation of microarray results (n = 5). All differences 
indicated between fasted and control fed are significant with a 
p-value< 0.05 (two-tailed, unpaired student's t-test). (B) Scatterplot 
of expression levels measured with qPCR and microarrays. Goodness 
of fit (r 2 ) is calculated with linear regression. Grey shading is according 
to legend in A. 



as functional p53 target gene [78]. In all three tissues 
investigated, Ddit4 mRNA is upregulated at latest by 
24 hours after onset of fasting and overnight fasting is 
sufficient to increase Ddit4 protein levels (Figure 6A), 
which has been shown by others in gastrocnemius 
muscle of rats [79]. In our data, differences in the magni- 
tude of fasting-induction between the mRNA and protein 
level (compare SM mRNA at the 24 h time point with 
protein levels in Figure 6A) could be explained by the fact 
that Ddit4 protein stability is highly regulated in different 
cell systems [80-82]. However, especially in adipose tissue 
Ddit4 protein seems to be stably induced when mice are 
fasted. To show that Ddit4 can be directly induced by p53 
in adipocytes, we treated mature C3H10T1/2 adipocytes 
for 6 hours with Nutlin-3, a specific p53 activator [83]. 
Nutlin-3 treatment led to an increase of Ddit4 mRNA 
similar to the canonical target Cdknla (Figure 6B). Most 
importantly, Ddit4 protein was also stably induced in all 



replicates (Figure 6B). In addition, the p53 targets Sesn2 
and Srebfl were regulated by Nutlin-3 in a way similar 
to the in vivo fasting situation (compare their expres- 
sion in Figure 6B and 4D). Hence, Ddit4 is stably in- 
duced by fasting and upregulated by p53 activation in 
cultured adipocytes. 

Overexpression of Ddit4 is sufficient to increase lipolysis 
in cultured adipocytes 

In a recent report Ddit4 was shown to be involved in lipid 
metabolism in adipocytes signaling via the mTORCl path- 
way [81]. Also in other studies, Ddit4 has been repeatedly 
described as a negative regulator of mTORCl in a variety 
of cell types [84,85]. Interestingly, in the context of 
starvation, the nutrient-sensitive mTORCl pathway 
needs to be suppressed for the proper fasting response 
in liver [86] and its suppression induces lipolysis in ad- 
ipocytes [87]. Hence, we examined whether upregula- 
tion of Ddit4 promotes lipolysis in adipocytes by 
inhibiting mTORCl activity. For this we transiently 
overexpressed Ddit4 in differentiated C3H10T1/2 adi- 
pocytes (Figure 7A) and determined glycerol and FFA 
in the medium as a measure of lipolysis. Indeed, we 
observed a -30% increased glycerol release and a -40% 
increased FFA release from Ddit4 overexpressing cells 
compared to the empty vector control (Figure 7B), 
while expression of genes in the lipolytic pathway 
remained unchanged (Figure 7A). Upon (3-adrenergic 
stimulation (1 hour 1 uM isoproterenol) the increase 
in lipolysis upon Ddit4 overexpression was still evi- 
dent, although not statistically significant (Figure 7C). 
However, as assayed by phosphorylation of the down- 
stream target S6K1 at threonine 389, mTORCl activity 
was unchanged despite effective overexpression of Ddit4 
(western blot in Figure 7A). An antibody against total 
S6K1 protein served as loading control and cells treated 
with rapamycin, a potent exogenous mTORCl inhibitor, 
as control for phosphorylation-specific S6K1 antibody. 
Others have reported that, in SM, dexamethasone- 
mediated Ddit4 increase leads to reduced mTORCl sig- 
naling [88], but, judging from our data, in adipocytes 
Ddit4-mediated lipolysis seems to be independent of 
mTORCl activity. This is consistent with the observation 
that the Ddit4-mTORCl axis is functional in some cell 
types but not in others [75]. Further, as knock-down of 
Ddit4 was reported to decrease insulin-stimulated de- 
novo lipogenesis in adipocytes [81], we wanted to rule out 
that the observed increase in lipolysis upon Ddit4 overex- 
pression is merely an effect of increased lipogenesis and 
with that higher lipid content per se. In addition to an 
unchanged phenotype as shown by mRNA expression of 
adipocyte-specific genes (Figure 7A), we did not detect an 
increase in the incorporation of radio-labeled glucose into 
total lipids in Ddit4 overexpressing cells compared to 
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Figure 6 Ddit4 is induced during fasting and upregulated by p53 activation in cultured adipocytes. (A) In vivo regulation of Ddit4 in 
fasted and fed mice. Relative expression values from qPCR measurements (n = 5) were normalized to 36b4 expression and the value for the first 
time point was set to 1 . Significance was determined with a 2-way ANOVA followed by a Bonferroni posttest to determine significance for the single 
time points (*p < 0.05; **p < 0.01; ***p < 0.001). Western blots were performed with tissue samples from six month old mice fed a chow diet 
(fed) or fasted overnight (fasted). Samples are from three littermates each. The following proteins served as loading controls (L.C.): (3-actin for 
WAT and LIV, (3-tubulin for SM. (B) 10 uM Nutlin-3, a specific p53 activator, was used to treat mature C3H10T1/2 adipocytes (day 7 of differentiation) for 
6 hours. qPCR mRNA measurements (n = 3) were normalized to 36b4 expression and related to control cells (treated with DMSO). A two-tailed, 
unpaired student's t-test was used to determine statistical significance (*p < 0.05; **p < 0.01; ***p < 0.001). The western blot shows data from 
three independent replicates treated with Nutlin-3 (+) or DMSO (-) using (3-actin as loading control. 



control cells (Figure 7D). Rather, a small decrease in lipo- 
genesis was observed (Figure 7D), which is consistent with 
a potential role of Ddit4 in regulating fasting responses in 
adipocytes. Hence, our data on the p53-target Ddit4 pre- 
sents a valuable example for hypothesis generation from a 
large-scale data set by suggesting a new role in the fine- 
tuning of the fasting-response in adipose tissue. However, 
further investigation is needed to elucidate the exact 
mechanism by which Ddit4 is enhancing lipolysis in adi- 
pocytes, because we show that this is independent form 
mTORCl and de-novo lipogenesis. 

Conclusions 

In this study, we took a comprehensive view on the fast- 
ing process in mice. Our combination of focused and 
genome-wide approaches reveals several fasting-related 
findings: (i) We provide a novel view on the immediate 
and dynamic response to fasting in mice during a 



48 hour period. These experiments focus on the timely 
regulation of liver genes relayed through the Ppara sig- 
naling pathway, which manifests in coordinated changes 
of serum parameters. The observed responses occur 
very early (between 3 to 6 hours) after onset of fasting 
and show simultaneous activation of different pathways, 
(ii) To our knowledge this is the first study focusing on 
the transcriptome response of white adipose tissue in 
fasted mice. With our bioinformatic analyses we identify 
an upregulation of apoptosis-related transcripts as well 
as a strong enrichment of transcriptional control com- 
ponents in the set of upregulated genes, (iii) Concen- 
trating our analyses on genes regulated in the three 
tissues mainly responsible for energy homeostasis during 
fasting (i.e. white adipose tissue, liver, skeletal muscle), the 
p53 signaling pathway appears to be a common and cen- 
tral regulator of fasting, possibly partly mediating its effect 
by down-regulation of the Srepb pathway, (iv) Finally, we 
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Figure 7 Transient overexpression of Ddit4 is sufficient to induce lipolysis in cultured adipocytes. (A)-(D) C3H10T1/2 cells were 
differentiated to adipocytes for 7 days. Cells were detached and electroporated in the presence of either empty vector (pMSCV) or Ddit4 
overexpression vector (pDdit4). (A) 48 hours later cells were harvested to measure mRNA (left qPCR) and protein (right western blots), and 
supernatants were collected for glycerol and FFA assay. qPCR measurements (n = 3) were normalized to 36b4 expression and related to control 
cells and a two-tailed, unpaired student's t-test was used to determine statistical significance (***p < 0.001). Western blot shows 2 biological 
replicates and a control treated with 100 nM rapamycin (Rapa). (B) and (C) Glycerol and FFA in cell supernatants were determined, normalized to 
protein content and related to measurements from control cells to yield the relative release caused by Ddit4 overexpression in unstimulated 
(B) and isoproterenol-stimulated (1 uM, 1 hour) (C) conditions. Significance was determined by a one-sample t-test (* p < 0.05). (D) Electroporated 
C3H10T1/2 adipocytes were incubated overnight with medium containing 14 C-labeled deoxy-glucose. Incorporated radioactivity was counted in 
the lipid extract and related to controls (n = 3). A statistically significant difference was not detected by a one-sample t-test. 



performed experiments that prove that Ddit4, a p53 target 
gene upregulated by fasting in all three tissues, is induced 
by p53 activation and sufficient to increase lipolysis in cul- 
tured adipocytes. In conclusion, our transcriptome study 
of three tissues combined with bioinformatic analyses and 
mechanistic in vitro experiments, suggests the p53-Ddit4 
axis as a novel mechanism in the fine-tuning of fasting 
common to major metabolic tissues. 

Methods 

Mouse experiments 

Experimental animal procedures were in accordance with in- 
stitutional guidelines and regulations of the University of 
Pennsylvania. The institutional review board of the Univer- 
sity of Pennsylvania reviewed and approved all mouse experi- 
ments. Male wild-type C57B1/6 J mice (Jackson Laboratories) 
were kept on regular chow (Research Diets) under standard 
conditions in a 12 hours day/ 12 hours night cycle. At an age 
of 10-12 weeks, animals were separated into two groups of 
25 mice each. Food was withdrawn from the fasting group at 



9 a.m., while the control group had continuous ad libitum 
access to their diet. Blood glucose was always determined 
before sacrificing mice. Mice were sacrificed at the beginning 
of the study (0 hours) and at 3, 6, 12, 24, and 48 hours after 
food removal. Mice were sacrificed by C0 2 inhalation and 
blood was collected by intra-cardiac puncture, serum isolated 
and stored at -80°C. Liver, epididymal adipose tissue, and 
skeletal muscle of the thigh were dissected in that order, flash 
frozen in liquid N 2 and stored at -80°C until mRNA extrac- 
tion. For western blot of Ddit4 tissues were dissected from 
six months old, male C57B1/6 J mice that were fasted 
overnight or continuously kept on normal chow diet. 

Serum parameters 

Blood glucose was measured by standard glucose oxidase 
glucometer test strips (One Touch Ultra, Fisher). Serum 
samples were analyzed using commercially available kits 
for insulin (#90080 Ultra Sensitive Mouse Insulin ELISA 
Kit, Crystal Chem), NEFAs (#995-34791 and 993-35191, 
Wako), glycerol (#F6428, Sigma), and p-hydroxybutyrate 
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(#2440-058, Stanbio). Corticosterone levels were de- 
termined with a Mouse/Rat Corticosterone ELISA kit 
(#55-CORMS-E01, Alpco). 

Tissue isolation 

Tissues were homogenized using a TissueLyser (Qiagen). 
mRNA was isolated by RNeasy spin columns and the 
RNeasy Lipid tissue kit (#74804, Qiagen), if needed. For 
tissue western blot, tissues were homogenized (Ultra- 
turax, IKA) in RIPA buffer, incubated 20 min on ice, and 
centrifuged (1000 g, 10 min, 4°C). Protein phase was iso- 
lated and measured with BCA kit (#23227, Pierce). 

qPCR analyses 

For qPCR measurements of tissue gene expression, iso- 
lated total RNA was reverse transcribed using High 
Capacity RNA-to-cDNA Master Mix (#4377474) and 
amplified using TaqMan Universal PCR Master Mix 
(#4324018) and measured using gene-specific Assays- 
on-demand (all Applied Biosystems). Amplifications 
were performed on an ABI Prism 7900HT machine fol- 
lowing manufacturers protocols. PCR efficiency was 
calculated from standard curves and the expression of 
36b4 (time series measurements) or Gapdh (verification 
of microarray measurements, Figure 5) was used for 
normalization. For cell culture experiments SYBR green 
qPCR was used. Total RNA was isolated with the 
GeneElute Mammalian Total RNA kit (#RTN70, Sigma). 
For reverse transcription Qiagen QuantiTect RT kit 
(#205311) was used. cDNA was then amplified using 
Sybr QPCR supermix (#11733-038, Life Technologies) on 
an ABI 7000 sequence detection system. Primers used 
Ddit4 (Fw: CCTGCGCGTTTGCTCATGCC; Rev: GGCC 
GCACGGCTCACTGTAT); Cdknla (Fw: GTCTGAGCG 
GCCTGAAGATTC; Rev: TGTTCCGGGCCCACCCGG 
GG); Sesn2 (Fw: CGCCACTCAGAGAAGGTTCA; Rev: 
ACGGGGTAGTCAGGTCATGT); Lpinl (Fw: GTCGT 
CGAGCAAGACAGATTCC; Rev: ACCAGGATCCCCA 
TTCTTGG); Srebfl (Fw: AAGCAAATCACTGAAGGAC 
CTGG; Rev: AAAGACAAGGGGCTACTCTGGGAG); 
Pparg2 (Fw: TGCCTATGAGCACTTCACAAGAAAT; 
Rev: CGAAGTTGGTGGGCCAGAA); ATGL (Fw: GTC 
CTTCACCATCCGCTTGTT; Rev: CTCTTGGCCCTCA 
TCACCAG); HSL (Fw: CCATCTCACCTCCCTTGG; Rev: 
TCCTTCCCGTAGGTCATAGG); Plinl (Fw: GGTACAC 
TATGTGCCGCTTCC; Rev: CTTTGCGCTCCGCCTCT); 
36b4 was used for normalizaiton (Fw: CAACCCAGCT 
CTGGAGAAAC; Rev: CCAACAGCATATCCCGAATC). 
Expression values were calculated with an in-house tool 
[89] employing the AnalyzerMiner algorithm [90]. 

Microarray experiments 

Frozen samples of white adipose tissue, liver, and skeletal 
muscle collected at the 24 hour time point were used for 



microarray experiments. RNA samples were quantified 
using a NanoDrop (ND-1000). RNA integrity was exam- 
ined using an Agilent 2100 Bioanalyzer. RNA samples 
(150 ng) with RNA integrity number >7 were used for 
target amplification and labeling via the Ambion WT 
Expression kit (#4411974) and Affymetrix WT Terminal 
Labeling kit (#900671) following manufacturers protocol. 
Mouse Gene 1.1 ST Array Plates (#901418, Affymetrix) 
were used for microarray hybridization, wash, stain and 
scan with GeneTitan hyb-wash-stain kits (#901622, 
Affymetrix) and a GeneTitan instrument. GeneTitan 
scanner data were collected with default parameters and 
further analyzed using Partek Genomics Suite. Data 
were normalized using default RMA method. A two-way 
ANOVA model with an interaction term between diet (fed 
or fasted) and tissue (liver, skeletal muscle and white adi- 
pose tissue) was set up. Pairwise comparisons were made 
between fed and fasted diet for each tissue. The resulting 
p-values of significance were corrected for multiple testing 
using Benjamini-Hochbergs false discovery rate (FDR) 
method. Genes within 5% FDR and changed at least by 
1.3-fold in either direction were called differentially 
expressed. Data was deposited in NCBI gene expression 
omnibus (GEO) with the accession number GSE46495. 

Functional annotations and mappings 

For DAVID functional annotation, Gene IDs of differen- 
tially regulated gene lists (1.3x, FDR5) were submitted to 
the DAVID website [43,91]. GO_FAT terms and KEGG 
pathways were considered significantly enriched if the 
Benjamini-Hochberg corrected p-value was >0.05. For 
functional clustering only enrichment scores (negative loga- 
rithm of geometric mean of p-values of entities in each 
cluster) larger than three were considered because they 
contain at least one entity with a significant p-value after 
multiple testing correction (<0.05). Gene-set enrichment 
analysis [45] was performed with all genes from the WAT 
microarrays as "expression data set" and the lists of tran- 
scription factor sub-classes, cofactors, and chromatin remo- 
delers from AnimalTFDB [44] as "gene sets database". 
The Venn diagram for intersection of genes differentially 
expressed in all tissues was drawn using Bio Venn [92]. 
MetaCore enrichment analysis matches gene IDs of 
possible targets for the "common", "similar" and "unique" 
sets (differentially expressed in all 3, 2 out of 3, and only 1 
tissue, respectively) with gene IDs in functional ontologies 
in MetaCore. The probability of a random intersection be- 
tween a set of IDs the size of target list with ontology en- 
tities is estimated in p-value of hypergeometric intersection. 
Heatmaps for Figure 4 were generated with Genesis [93]. 

cDNA cloning 

The Ddit4 coding sequence was PCR-amplified with cod- 
ing sequence-flanking primers from mouse adipose tissue 
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cDNA and cloned into a pMSCV mammalian expression 
vector (Life Technologies) between Xhol and EcoRI re- 
striction sites using standard procedures. Correct cloning 
was verified by sequencing of the whole insert. 

Cell culture experiments 

C3H10T1/2 cells were maintained in growth medium 
(High-glucose Dulbeccos modified Eagle s medium sup- 
plemented with 10% FBS, 2 mM L-glutamine, 100 U/ml 
penicillin, 100 mg/ml streptomycin (all from Life Tech- 
nologies)). Two days post-confluent cells were induced to 
undergo adipogenesis by addition of 1 \iM dexametha- 
sone, 500 [iM 3-isobutyl-l-methylxanthine, 5 (ig/ml 
insulin (all Sigma), and 1 \iM rosiglitazone (Alexis) as 
described by others [94]. From day 3 on growth 
medium was only supplemented with 1 (ig/ml insulin 
for 2 days before switching back to normal growth 
medium. For Nutlin-3 treatments day 7 adipocytes were 
treated for 6 hours with 10 \xM of Nutlin-3 (Sigma) or 
DMSO as control, before cells were harvested for RNA 
and protein analysis. For Ddit4 overexpression day 7 adi- 
pocytes were detached with a trypsin (0.25%)/collagenase 
(0.5 mg/ml) mix, washed and resuspended in electropor- 
ation buffer R (Neon electroporation kit, #MPK1025, Life 
Technologies) containing 1 \ig of either empty overexpres- 
sion vector (pMSCV, Life Technologies) or vector with 
Ddit4 coding sequence at a concentration of 30,000 cells/ ul. 
Electorporation was performed in 10 \A tips with 1400 V/ 
30 ms pulses using a Neon transfection system (Life Tech- 
nologies). Three electroporation reactions were reseeded 
in one well of a 12-well plate in growth medium without 
antibiotics, which was replaced by normal growth medium 
on the next day for RNA, protein and glycerol measure- 
ments. For FFA measurements 2% FFA-free BSA (#K31- 
002, PAA) containing growth medium without FBS was 
used. These media were supplemented with 100 nM rapa- 
mycin or 1 \iM isoproterenol for 1 hour where indicated. 
After 48 hours, cells were harvested for protein and RNA 
analysis, and supernatants were collected for glycerol and 
FFA determination. 

Free fatty acid and glycerol measurements 

Collected media of electroporated cells were centrifuged 
at 12,000 g for 5 min and supernatants were transferred 
to new tubes. FFA and glycerol contents were measured 
using commercial kits (from Wako (#NEFA-HR(2)) and 
Thermo Scientific (#TR22421), respectively) according 
to the manufacturer s protocol. Concentrations were de- 
rived from standard values/curves and related to the 
amount of protein in the same well as determined with 
BCA assay (#23227, Pierce). Finally, values were related 
to the empty vector control measurements to obtain 
relative glycerol/FFA release in per cent. 



Lipogenesis assay 

C3H10T1/2 adipocytes were incubated overnight with 
medium (no-glucose DMEM, 10% FBS, 50 ug/mL 
streptomycin, 50 units/mL penicillin) supplemented 
with 0.5 g/L glucose and 0.1 uCi D[ 14 C(U)] -glucose/ml 
(ARC). Cells were washed three times with ice-cold 
PBS before cellular lipids were extracted with hexane/ 
isopropanol (3:2, vol). The incorporated radioactivity 
in the organic phase was determined by liquid scintilla- 
tion counting. Counted values were corrected by pro- 
tein content. 

Western blot analysis 

Western blot analysis was performed as we described 
previously [95] with the following changes: SDS -lysis 
buffer was supplemented with phosphatase inhibitor 
(PhosStop, Roche) to reduce changes in phosphorylation 
states, a 4-12% Bis-Tris gel (Life Technologies) was used 
and 40 \ig protein was loaded per lane. Antibodies used: 
a-Ddit4 (# 10638-1-AP, ProteinTech Europe), a-S6Kl 
(#9202, Cell Signaling), a-phospho S6K1 (Thr389) (#9206, 
Cell Signaling). For tissue western blot 70 \ig of protein 
was loaded. The following antibodies were used to detect 
loading controls: p-actin (#A5316, Sigma) for WAT and 
LIV and (3-tubulin (#T5201, Sigma) for SM. Detection 
was performed using ECL prime substrate from GE 
Healthcare. Before reprobing blots were stripped with 
Restore WB stripping buffer from Pierce. 

Statistical analyses 

To determine statistical significance in time series mea- 
surements a 2-way ANOVA was used followed by a 
Bonferroni posttest to determine significant differences 
for the single time points. For comparative qPCR mea- 
surements upon Ddit4 overexpression and Nutlin-3 treat- 
ment a two-tailed, unpaired students t-test was used and 
for FFA, glycerol and lipogenesis measurements a one- 
sample t-test. A p < 0.05 was considered as statistically 
significant (*p < 0.05; **p < 0.01; ***p < 0.001). qPCR 
time series measurements and bar graphs are shown as 
average ± SEM from independent experiments (sample 
size as indicated in figure legends). Significance compu- 
tation in the microarray data was performed using a 
two-way ANOVA model with a Benjamini-Hochbergs 
false discovery rate (FDR) to correct for multiple testing. 
In this study, transcripts with an FDR5 and > 1.3-fold 
difference were considered as differentially expressed be- 
tween fasted and fed groups. For DAVID analyses GO 
terms and KEGG pathways were considered as signifi- 
cantly enriched if the Benjamini-Hochbergs corrected p- 
value was <0.05 [43]. Significant enrichment of gene lists 
in transcription factors, cofactors, and chromatin remode- 
lers was tested with a Chi-square test with Yate s correc- 
tion (http://graphpad.com/quickcalcs/contingencyl.cfm). 
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Additional files 



Additional file 1: Lists of probes expressed more than 1.3-fold with 
an FDR5 for each of the three tissues; p-value given is corrected for 
multiple testing. 

Additional file 2: Detailed results of the functional clustering by 
DAVID including terms/pathways contributing to each cluster; Lists 
of genes up- and downregulated by 24 hours of fasting in WAT 
where submitted separately and resulting clusters are listed in 
separate tabs. 

Additional file 3: Annotation of genes regulated in fasting WAT 
according to AnimalTFDB; Results of GSEA analysis testing for 
enrichment in AnimalTFDB categories (sub-categories of 
transcription factors (TF), cofactor (CoF), chromatin remodelers 
(ChrRem)); Manual mapping of differentially expressed WAT genes 
in the AnimalTFDB categories TF, CoF, and ChrRem; GSEA results 
for gene sets enriched in up- or downregulated genes (by fasting in 
WAT) using c2 and c5 (MSigDB) modules as gene sets. 

Additional file 4: Metacore analysis focusing on networks 
overrepresented in the common list reveals the p53 node as the 
second highest scoring network hub. 

Additional file 5: Metacore pathway analysis on Affymetrix-probes 
regulated commonly by fasting in WAT, LIV, and SM yields 
Srebp-regulated pathways as 2 top-scoring entities. 

Additional file 6: DAVID mapping to GO_BP_FAT of up- or down- 
regulated genes for WAT, LIV, and SM; only significant terms are 
depicted. 
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